Bimodal behavior of the heaviest fragment distribution in projectile fragmentation 
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The charge distribution of the heaviest fragment detected in the decay of quasi-projectiles pro- 
duced in intermediate energy heavy-ion collisions has been observed to be bimodal. This feature is 
expected as a generic signal of phase transition in non-extensive systems. In this paper we present 
new analyses of experimental data from Au on Au collisions at 60, 80 and 100 MeV/nucleon show- 
ing that bimodality is largely independent of the data selection procedure, and of entrance channel 
effects. An estimate of the latent heat of the transition is extracted. 

PACS numbers: 05.70.Fh Phase transitions: general studies ; 25.70.-z Low and intermediate energy heavy- 
ion reactions ; 25.70.Pq Multifragment emission and correlations 



At a first-order phase transition, the distribution of the 
order parameter in a finite system presents a character- 
istic bimodal behavior in the canonical or grandcanon- 
ical ensemble 0, S 0, II- The bimodality comes from 
an anomalous convexity of the underlying microcanonical 
entropy 0. It physically corresponds to the simultane- 
ous presence of two different classes of physical states for 
the same value of the control parameter, and can survive 
at the thermodynamic limit in a large class of physical 
systems subject to long-range interactions In the case 
of nuclear multifragmentation, a natural order parame- 
ter is the size of the heaviest cluster produced in each 
collision event. Indeed this observable provides an order 
parameter for a large class of transitions or critical phe- 
nomena involving complex clusters, from percolation to 
gelation, from nucleation to vaporization, from reversible 
to irreversible aggregation [3, 0, S] • 

In this context, the recent observation by the INDRA- 
ALADIN collaboration [l(| of a sudden change in the 
fragmentation pattern of Au quasi-projectiles, loosely re- 
ferred to as bimodality, has triggered a great interest in 
the heavy- ion community TJ|. Looking at the correla- 



tion between the two heaviest fragments emitted in each 
event as a function of the violence of the collision, a clear 
transition is observed between a dominant evaporation- 
like decay mode, with the biggest cluster much heav- 
ier than the second one, and a dominant fragmentation 
mode, with the two heaviest fragments of similar size. 
A similar behavior has been reported in ref. Q. Differ- 
ent physical scenarios have been invoked to interpret the 
phenomenon: finite-system counterpart of the nuclear 
matter liquid-gas phase transition lid. 



transition of highly deformed systems TJ] , self-organized 
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criticality induced by nucleon-nucleon collisions T^, It 
In [9|], the two decay modes were associated to different 
excitation energies, suggesting a temperature-induced 
transition with non-zero latent heat. The qualitative 
agreement between refs. 0, [l(| suggests that bimodality 
is a generic phenomenon. However, differences between 
the two data sets subsist, and trigger or selection bias 
cannot be excluded. To disentangle between the differ- 
ent scenarios, it is necessary to control the role of the 
entrance channel dynamics and establish if the transi- 
tion is of thermal character. In this letter, to progress 
on these issues, event ensembles with equiprobable exci- 
tation energy distribution are built and compared. 

We present a new analysis of quasi-projectiles (QP) 
produced in Au+Au collisions measured with the IN- 
DRA apparatus [T3| at the GSI laboratory at incident 
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energies from 60 to 100 MeV/nucleon [l8(. The robust- 
ness of the signal of bimodality is tested against two dif- 
ferent QP selection methods. A weighting procedure [l3j 
is applied to test the independence of the decay from 
the dynamics of the entrance channel. Finally, a double 
saddle-point approximation is applied to extract from the 
measured data an equivalent-canonical distribution that 
can be quantitatively confronted to statistical theories of 
nuclear decay flp| ] - 

In this energy regime, a part of the cross section corre- 
sponds to collisions with dynamical neck formation 21 1. 
We thus need to make sure that the observed change 
in the fragmentation pattern jTo[ is not trivially due to a 
change in the size of the QP. After a shape analysis in the 
center of mass frame 20] , only events with a total forward 
detected charge larger than 80% of the Au charge were 
considered (quasi-complete events). Two different pro- 
cedures aiming at selecting events with negligible neck 
contribution were adopted. In the first one [10( (I) by 
eliminating events where the entrance channel dynamics 
induces a forward emission, in the quasi-projectile frame, 
of the heaviest fragment Z\ [22| ■ For isotropically decay- 
ing QPs, this procedure does not bias the event sam- 
ple but only reduces the statistics. In a second strategy 
(II) the reduction of the neck contribution is obtained by 
keeping only "compact" events by imposing (i) an upper 
limit on the relative velocity among frag ments, and (ii) 
a QP size constant within 10%, see |12|_for details. In 
both cases fission events were removed [lOj. 

The selected samples contain altogether about 30% of 
the quasi-complete events at the three bombarding ener- 
gies. The main characteristics of the distribution of the 
heaviest fragment are presented in Fig. [TJ as a function 
of the total transverse energy of light-charged products 
[Z = 1,2) 12311. An excitation energy scale, estimated by 
calorimetry [24j 28], is also given. 

For increasing violence of the collision, the average size 
of the largest fragment monotonically decreases. The av- 
erage behavior is smooth, but higher moments of the dis- 
tribution reveal a clear change from the high Z\ evapo- 
ration dominated pattern, to the low Z\ multifragmen- 
tation dominated one, passing through a region of max- 
imal fluctuations where the skewness changes its sign. 
These moments appear relatively independent of the se- 
lection criterion. About one event out of four is common 
between the two sets; the differences in the observables 
evaluated with the two criteria thus give an estimation 
of the bias induced by the selection of data. The rela- 
tive abundances observed in the correlation between the 
charge of the heaviest fragment and the deposited excita- 
tion energy are clearly governed by the impact parame- 
ter. The presence of a sudden jump in the most probable 
Z\ value depends on the selection method and cannot be 
taken as a signature of a transition, as it was proposed in 
previous works 2, 14, 15, lij]. The only veritable proof 
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FIG. 1: Upper part: average (dots), standard deviation 
(squares) and skewness (triangles - right Y-axis) of the dis- 
tribution of the heaviest fragment as a function of the light- 
charged particles transverse energy at an incident energy of 
80 MeV/nucleon. Lower part: correlation between the charge 
of the heaviest fragment and the calorimetric excitation en- 
ergy. The open squares indicate the most probable Zl values. 
The average total source size Z s is given by the full line. Left 
side: selection (I); right side: selection (II). 



bumps in the Z\ distribution for a system in thermal 
contact with a heat reservoir at the transition tempera- 
ture [J, 0, H, IH • However, the distribution of the energy 
deposit in a heavy-ion collision is not determined by ran- 
dom exchanges with a thermal bath. This means that the 
experimental ensemble is not canonical and the Z\ distri- 
bution has no meaning in terms of statistical mechanics. 
To cope with this problem, a simple procedure has been 
proposed in ref. [13j. The bimodality in the canonical 
two-dimensional probability distribution pp(E* , Z\) of a 
system of given size Z a at a first order phase transition 
point reflects the convexity anomaly of the underlying 
density of states W Zs {E*,Z{) [H, i, % according to: 



pp(E*, Zi) - W Zs {E\ Z x ) ex V (-pE*)Z £ 



(1) 



where Zp is the partition function. In an experimental 
sample, the energy distribution is not controlled by an 
external bath through a Boltzmann factor, but it is given 
by a collision and detector dependent functional g(E*): 



p exp (E*)<x j dZ 1 W Zs (E*,Z 1 )g(E*). 



(2) 



of bimodality would be the observation of two distinct 



The convexity of the density of states can be directly 
inferred from the measured experimental distribution, by 
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a simple weighting of the probabilities associated to each 
deposited energy: 



p w {E*,Z l ) = 



p exp {E*,Z x ) _p f3 (E*,Z 1 ) _ W Ze {E*,Zx) 



Pexp(E*) pp{E*) Wz B {E*) 

(3) 

This procedure allows to get rid of the entrance chan- 
nel impact parameter geometry that naturally favors the 
lower part of the E* distribution. To produce a flat E* 
distribution according to eq.([3]), we have weighted the Z\ 
yields in each E* bin with a factor proportional to the 
inverse of the bin statistics. 

The results obtained with the two different selection 
methods are given in Fig. [5] (bottom) . To take into ac- 
count the small variations of the source size, the charge 
of the heaviest fragment Z\ has been normalized to the 
source size. After the weighting procedure, a bimodal 
behavior of the largest fragment charge clearly emerges 
in both cases. 
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FIG. 2: Upper part: measured distribution of the charge of 
the largest fragment normalized to the charge of the source 
detected in Au+Au collisions at three different bombarding 
energies. Lower part: normalized distributions obtained con- 
sidering the same statistics for each excitation energy bin. 
The left (right) side shows distributions obtained with the 
data selection method (I) ((II)). 



Eq. j3|) holds only if the bias function g in eq. ((2|) does 
not explicitly depend on Z\ , which is a phase-space domi- 
nance assumption. The physical meaning of this hypoth- 
esis is that the entrance channel geometry and dynam- 
ics (as well as the bias induced by the detection system 
and data selection) determine only the energy distribu- 
tion and size (Z s ) of the QP, while for each given value 
of E* and Z s the size of the heaviest fragment (Zi) is 



dominated by the corresponding available phase space. 
The similarity of the two samples at 80 MeV/nucleon, 
after the weighting procedure, is an indication that the 
bias induced by the data sorting is small. The phase- 
space dominance hypothesis can further be checked by 
comparing the effect of the weighting procedure on data 
issued from different entrance channel dynamics. This 
is done in Fig. [2j where the same weighting method has 
been applied on data at different bombarding energies. 
The comparison is not conclusive in the case of selection 
(I) , where the excitation energy distributions obtained at 
the different incident energies happen to be largely super- 
posable (Fig. [2] top left), and we cannot a priori exclude 
a bias function. Conversely in the case of selection (11), 
we can see that the weight of the low Z\ component, as- 
sociated to more fragmented configurations and higher 
deposited energy, increases with the bombarding energy. 
This difference disappears when data are weighted, show- 
ing the validity of the phase-space dominance hypothesis. 

The three studied energies and the two selection crite- 
ria (I) and (11) produce similar but not identical distribu- 
tions even after renormalization, meaning that a residual 
bias on the density of states exists. One may ask whether 
this bias prevents a sorting and dynamic-independent ex- 
traction of the entropic properties of the system. To an- 
swer this question, we can compare the information on 
the coexistence zone in the {Z\, E*) plane extracted from 
the different samples. We thus have to solve eq.Q for 
the canonical distribution pp t (E* , Z\) at the transition 
temperature (3 t at which the two peaks of the energy dis- 
tribution have the same height [4[ . This is easily obtained 
in a double saddle point approximation [l3j |: 



E 

i=l, g 



N, 



1 



V detYi. 



exp 



1 



(4) 



where xi = (E* — E\, Z\ — Zi). £j represents the 
variance-covariance matrix and is related to the entropy 
curvature matrix (see formulae 10, 11 and 12 of (H}). The 
correlation coefficient p=oziE' / o~ZiO~e* , which is one of 
the parameters, was calculated from the data at the three 
incident energies, before the weighting procedure and for 
each selection method, on the largest validity domain i.e. 
1-8 MeV/nucleon for (I) and 1-12 MeV/nucleon for (II) 
(see tabic fl}. Ej is evaluated at the liquid I (gas g) solu- 
tion, and Ni are the proportions of the two phases, with 
Ni/N g = VoMTTi/y/MYr g . 

The weighted experimental distribution can be fit- 
ted with the function p w (E*, Z\) = pp t (E*, Z{)/pp t (E*) 
which, using eq.([4]), is an analytic function, p being 
fixed, we have performed an 8-parameter fit with the 
two data sets corresponding to the two selection proce- 
dures at the two higher bombarding energies on the ex- 
citation energy range 2-7 MeV/nucleon; to avoid small 
number effects only 2D-bins with significant statistics 
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TABLE I: Parameters of the equivalent canonical distribution eq.(|4]) at the transition temperature as estimated from the two 
data selection methods. The \ 2 °f the fit is also given. 
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X 2 /N dof 


set (I) E/A= 


=80MeV 


-0.861 


72.5 


16.5 


1.42 


2.25 


12.1 


13.4 


8.52 


2.62 


0.53 


set (I) E/A= 


lOOMeV 


-0.861 


69.3 


15.9 


1.67 


2.30 


12.1 


13.7 


8.76 


2.83 


0.59 


set (II) E/A= 


=80MeV 


-0.925 


69.1 


12.6 


1.02 


1.78 


2.10 


24.6 


10.4 


4.04 


0.80 


set (II) E/A= 


=100MeV 


-0.925 


68.3 


12.5 


1.07 


1.77 


2.96 


24.4 


10.2 


3.96 


0.96 



(> 0.5% of the corresponding E* slice) were used. The 
obtained parameter values are given in table [U In par- 
ticular, we can estimate the latent heat of the transi- 
tion of the heavy nuclei produced as AE = E g — E{ = 
8.1(±0.4) stot .(+1.2- 0.9)ay St . MeV/nucleon. Statistical 
error was derived from statistical errors on Ei and E g 
and systematic errors from the comparison between se- 
lections (I) and (II). The latent heat is derived from a 
difference and so the possible effect of systematic errors 
in the determination of excitation energy by calorimetry 
due to detection limitations (neutrons are not detected 
nor fragment masses measured) (2{| should be included 
in given error bars. Note also that the deduced param- 
eter values Ei and E g are outside the excitation energy 
range used for the fit. 
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entrance channel effect of the impact parameter has been 
removed by weighting the Zi distribution by the statistics 
of the excitation energy distribution, a clear indication of 
bimodality in the decay pattern is observed. This behav- 
ior appears to be robust against the selection method, 
the entrance channel dynamics and the estimator of the 
deposited excitation energy. This analysis supports the 
interpretation of the discontinuity already observed in 
the decay pattern [l(| as the finite system counterpart 
of a first order phase transition. A multidimensional fit 
allows to extract, through a double saddle point approx- 
imation, the coexistence zone and a first estimate of the 
latent heat of the transition. 

The present results are coherent with other signals 
from Au quasi-projectiles considered indicative of a first 
order phase transition like a fossil signal of spinodal fluc- 
tuations and configurational energy fluctuations associ- 
ated with ne gat ive heat capacity [26 



271] . Interpretations 
given m do not register m that coherent pic- 

ture. However it would be interesting to know if those 
interpretations can verify the bimodality of Z\ for the 
weighted distribution and its independence of the inci- 
dent energy as it is observed in that work. 
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FIG. 3: Experimental distribution of the largest cluster charge 
normalized to have identical statistics for each excitation en- 
ergy bin, with the two different data selection techniques (I) 
and (II) and for 80 MeV/nucleon incident energy. Different 
estimators of the deposited excitation energies are considered. 

Finally we use other estimators such as the total for- 
ward charged product multiplicity M to t and the trans- 
verse energy Et\2- The measured distributions weighted 
via eq.© with these different estimators are presented 
in Fig. [3] We can see that bimodality is preserved in all 
cases, and the different energy estimators predict close 
positions for the two peaks. 

To conclude, in this paper we have presented a com- 
parative analysis of the quasi-projectile Au+Au data col- 
lected with the INDRA apparatus at incident energies 
between 60 and 100 MeV/nucleon. Two different meth- 
ods for quasi-projectile selection have been used, which 
do not select the same physical events. Once the trivial 
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